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ABSTRACT 

' The equations governing dissipative relativistic hydrodynamics are formu- 

lated within the 3+1 approach for arbitrary spacetimes. Dissipation is ac- 
counted for by applying the theory of extended causal thermodynamics (Israel- 



> 

Stewart theory). This description eliminates the causality violating infinite 



signal speeds present in the conventional Navier-Stokes equation. As an ex- 
| ample we treat the astrophysically relevant case of stationary and axisymmet- 

ric spacetimes, including the Kerr metric. The equations take a simpler form 

O 

Q" 1 whenever the inertia due to the dissipative contribution can be neglected. 

U 

bJ[), Key words: relativity - hydrodynamics - black hole physics - accretion, 



accretion discs - galaxies: active - stars: neutron 

1 INTRODUCTION 

The motion of dissipative fluids in strong gravitational fields is of considerable interest 
in various fields of astrophysics and cosmology. Examples include accretion discs around 
compact objects, rotating relativistic fluid configurations such as supermassive stars, neutron 
stars or strange (boson) stars, the collapse of stellar objects and the merging of compact 
objects. Examples in cosmology cover inflationary cosmological scenarios and the evolution 
of density fluctuations. 

The non-stationary modeling of relativistic matter is most conveniently performed within 
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2 J. Peitz and S. Appl 

the 3+1 formalism, where the equations of motion for gravitational- and matter fields are 
decomposed w.r.t. a congruence of fiducial observers (FIDOs), allowing to express time 
derivatives on a per-unit-universal-time basis. The 3+1 representation of the equations for 
ideal, non-dissipative (relativistic) hydrodynamics as well as Maxwell's equations for the 
electric and magnetic fields and Einstein's equations for the gravitational fields have been 
discussed by many authors, both in a cosmological context (e.g. Durrer & Straumann 1988) 
and in the case of black hole spacetimes (e.g. Thorne & Macdonald 1982). 

Modeling dissipative processes requires non-equilibrium or irreversible thermodynamics. 
Standard (or classical), irreversible thermodynamics (in the following referred to as standard 
thermodynamics) was first extended from Newtonian to relativistic fluids by Eckart (1940). 
However, the Eckart theory, and a variation thereof by Landau & Lifshitz (1959) shares with 
its Newtonian counterpart serious problems. Notably that dissipative fluctuations propagate 
at an infinite speed. In addition, generic short wavelength secular instabilities driven by 
dissipative processes exist (Lindblom & Hiscock 1983) and finally, no well-posed initial value 
problem exists for rotating fluid configurations. 

At the origin of these problems in standard thermodynamics is the description of non- 
equilibrium via the local equilibrium states alone, i.e. it is assumed that local thermodynamic 
equilibrium is established on an infinitely short time-scale (see e.g. Jou, Casas- Vazquez & 
Lebon 1997 for an introduction). In extended theories of irreversible thermodynamics the 
set of thermodynamic variables is extended to include the dissipative variables. This re- 
stores causality and stability under a wide range of conditions (Hiscock & Lindblom 1983). 
A non-relativistic extended theory was proposed by Miiller (1967), and was then general- 
ized to the relativistic case by Israel (1976) and Stewart (1977). The extended theory is 
commonly referred to as causal thermodynamics, second-order thermodynamics or transient 
thermodynamics. 

The problem of non-causality has recently received attention in the context of transonic 
accretion discs. Two different approaches have been proposed to overcome this difficulty. 
For steady flow Narayan (1992) has established causality by calculating the coefficient of 
kinematic viscosity within an extended version of flux limited diffusion theory (Levermore & 
Pomraning 1981), assuming a particular steady state phase-space distribution function for 
the turbulent fluid elements in the disc. The influence of this modified viscosity coefficient 
was studied in stationary accretion discs by Popham & Narayan (1994) and Syer & Narayan 
(1993). A relativitic generalization of the modified viscosity has been proposed, and used in 
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3+1 formulation of non-ideal hydrodynamics 3 



models of stationary relativistic accretion discs (Peitz & Appl (1997)). A different approach 
by Papaloizou & Szuskiewicz (1994), which is related to a causal description for the thermo- 
dynamics, was used by Kley & Papaloizou (1997) in time- dependent models for accretion disc 
boundary layers. Gammie & Popham (1997) have recently considered a similar extension to 
stationary relativistic accretion discs. In cosmology the theory of causal thermodynamics is 
currently attracting growing interest predominantly in the contexts of re-heating processes 
after inflation (Zimdahl, Pavon & Maartens 1997) and in linear perturbation theory for the 
evolution of density fluctuations (Maartens & Triginer 1997). 

This paper provides a complete set of equations for dissipative fluid mechanics in their 
3+1 representation, using a causal description of thermodynamics. In Section 2 the ba- 
sic elements of both standard and causal dissipative hydrodynamics are reviewed in their 
spacetime description, and their 3+1 representation is derived in Section 3. As a particular 
application of astrophysical interest we specify the system of equations given in Section 3 
to the case of a stationary, axisymmetric back ground in Section 4. 

2 DISSIPATIVE RELATIVISTIC HYDRODYNAMICS 

The equations of ideal relativistic fluid mechanics and the equations for dissipative relativis- 
tic fluid mechanics in both the standard irreversible and the extended causal thermodynam- 
ics description are reviewed. For a detailed discussion see Israel & Stewart (1979), Hiscock 
& Lindblom (1983) or a recent treatment by Maartens (1997). 

2.1 Notation 

We use geometrized units such that c — 1 — G. Tensor fields defined on spacetime (Ai,g) 
with metric g of signature (—,+,+,+) appear in roman (e.g. u, T), while scalar functions 
are in italic. The velocity u of the fluid is normalized to u • u = — 1. The tensor h = g + u <g> u 
projects into the 3-space orthogonal to u, the local rest frame of the fluid (LRF). Total 
projections (i.e. projection in any free index) parallel to and orthogonal to u are denoted 
by ( ) u and ( ) h . If A, B and C are a scalar, vector and rank-2 tensor field on (JA,g) , 
respectively, then 



A u = 



A h = A , 



B 



'u 



5, 



'u 



u-B 



B h = h • B 



C u = C u = u-C 



u , C h = h-C-h. 



(1) 
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4 J. Peitz and S. Appl 

Two covariant differential operators D and D are defined as projections of the affine con- 
nection V on (.M,g) into directions parallel and orthogonal to u, 

D = V u = u • V , (2) 
D = Vh = h • V . (3) 

For 2-tensor fields C we further introduce (anti-) symmetrization operators ( ) a and ( ) s by 
C = (C) a if C is anti-symmetric and C = (C) s if C is symmetric and trace-free. Finally, the 
irreducible decomposition of Vu yields 

Vu = cx+cu + l/36>h-a<g}u, (4) 

with shear cr = (Vu)^, vorticity cu = (Vu)^, acceleration a = Du and expansion = D-u = 
V • u of the fluid trajectories. 

2.2 Perfect fluids 

A perfect fluid is described by the velocity u, baryon number density n, mass-energy density 
p, isotropic pressure p and specific entropy s, which are subject to the conservation laws 

= V-n, (5) 
= V-T. (6) 

The particle current vector n and the symmetric stress-energy tensor T are given by 

n = nu , (7) 
T = pu <g> u + ph . (8) 

The LRF conservation laws for energy and momentum result from projecting (6) parallel 
and orthogonal to u. With (8) one can write = (V • T) u and = (V ■ T) h as 

= D p + ( p + p )e, (9) 
= Dp+(p + p)&. (10) 

The metric g is coupled to stress-energy T by Einstein's equations 

G = 8ttT, (11) 

where G is the Einstein tensor. 

Thermodynamic scalar functions are defined in the LRF. The entropy flux 

s = sn (12) 
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3+1 formulation of non-ideal hydrodynamics 5 
is conserved along flow lines (adiabatic flow), 

= V-s. (13) 
The temperature T is defined via the Gibbs equation 

Tds = d(p/n)+pd(l/n) , (14) 

where d is the exterior derivative on (Ai,g) . In general, two thermodynamic scalars are 
needed as independent variables, which we choose n and p. A scalar equation of state, e.g. 
p = pin, p), closes the system of equations (5), (9)-(ll) for the dynamical variables {n, p, u, 

g}- 

2.3 Dissipative fluids 

Choosing the particle current for dissipative fluids by (7) corresponds to selecting an average 
velocity u such that the particle flux in the associated rest frame vanishes. This so-called 
particle frame or Eckart frame (Eckart 1940) is the natural frame in systems where particle 
number is conserved (see Israel & Stewart 1979 for the alternative energy frame description). 
The state of the fluid is assumed close to a fictitious thermodynamic equilibrium state, 
characterized by the local thermodynamic equilibrium scalars n ,p ,p , so,T and the local 
equilibrium velocity Uo, which in the Eckart frame can be chosen such that only the pressure 
p deviates from the local equilibrium pressure po by the bulk viscous pressure II = p — p , 
whereas n = n and p = p . Dropping subscripts allows then to write the general stress- 
energy tensor for dissipative fluids as 

T = pu®u+ (p + i7)h + q®u + u(g)q + 7r, (15) 

where the heat flux q relative to the particle frame and the anisotropic stress tensor n are 
orthogonal to u, 

q = q h , n=(n) s h . (16) 

Conservation laws again hold for n and T. However, in irreversible thermodynamics the 
entropy is no longer conserved. According to the second law, the rate of entropy generation 
must therefore be positive definite, 

V-s > 0, (17) 

implying that s has a dissipative vector contribution R in excess to (12), 

s = sn + R/T. (18) 
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6 J. Peitz and S. Appl 

Following the phenomenological Israel-Stewart approach, s = Sq remains related to T = T 
by the Gibbs equation (14). The dissipative part R is assumed to be an algebraic function of 
n and T only, which vanishes in equilibrium (R = 0). The theories of standard irreversible 
thermodynamics and of extended causal thermodynamics differ in the forms of R, as given 
below. 

The equations of energy and momentum conservation for (15) can be written as 

= Dp + (p + p + n) 0+ a : n+ (2a + D) • q , (19) 
= (p + p + iI)a + D(p + iI) + D-7r + 7t-a+(i)q) h + (a+cu + 4/36>h)-q (20) 

where B : C = tr(B • C) for 2-tensor fields B, C. The above treatment applies for a single- 
component fluid and allows a natural extension to multi-component fluids. 

2.4 Standard thermodynamics 

Standard thermodynamics assumes a linear dependence of R on the thermodynamic fluxes 
{II, q, 7t}, which is possible only if (18) takes the form 

s = sn + q/T . (21) 
Using (5), (19) and (20) yields the entropy generation rate 

TV-s = -6>i7-(DlnT + a)-q-0:7T. (22) 

The simplest relation to be imposed between the thermodynamic fluxes {II, q, tc} and the 
thermodynamic forces {30, D InT + a, a} in agreement with (17) is linear, 

n = -co, (23) 

q = —AT (D InT + a) , (24) 
7t = -2r}U , (25) 

with non-negative coefficients of bulk viscosity ((p,n), thermal conductivity, \{p, n) and 
shear viscosity r)(p,n). This brings (22) to 

TV-s = i7 2 /C + q-q/(AT) +71:71/(277), (26) 

which on using (5), (14), (19) and (20) yields an evolution equation for the entropy density 

s, 

TnDs = -On- V-q-a-q-7t: a. (27) 
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3+1 formulation of non-ideal hydrodynamics 7 

As a consequence of (23)- (25), the thermodynamic fluxes react instantaneously to the cor- 
responding thermodynamic forces, implying propagation of signals at (causality violating) 
infinite speed. In standard thermodynamics thedynamics of the fluid is governed by the (com- 
pressible) Navier-Stokes equation, which results from substitution of (23)- (25) into (20). 

2.5 Causal thermodynamics 

Kinetic theory can motivate that R is second-order in the dissipative terms (Israel & Stewart 
1979). Truncation at first order removes terms necessary for causality and stability. The most 
general algebraic form for R of at most second-order in the dissipative fluxes leads to 

s = sn + q/T - 1/2 (f3 II 2 + ftq • q + fan : n) u/T + a IIq/T + a l7 r • q/T . (28) 

The entropy density measured in LRF then becomes 

-s-u = sn- 1/2 ((3 n 2 + (3 iq -q + (3 2 n: n) . (29) 

The negative sign of the non-equilibrium contributions reflects the fact that the entropy den- 
sity is maximum in equilibrium. The thermodynamic coefficients f3j(p,n) > in (28) model 
deviations of the physical entropy density from sn due to scalar /vector /tensor dissipative 
contributions to R. The aii(p,n) model contributions due to viscous/heat coupling, which 
do not influence the physical entropy density (29). 

The entropy generation rate associated with (28) follows from (5), (14), (19) and (20) to 



TV • s = -nX - q • Y - 71 : (Z) s (30) 
with the scalar, vector and rank-2 tensor fields 

X = Q + (3 L>n - a D • q - K Tq ■ V (a /T) + 1/277TV • ((3 u/T) , (31) 
Y = VlnT + a + / 9i J Dq-a Vi7-a 1 V-7t 

-(1 - K )i7TV (oq/T) - (1 - Kl )Tn ■ V K/T) + l/2TqV • (Au/T) , (32) 

Z = Vu + f3 2 Dn - e*i Vq - K t Tq ® V («i/T) + 1 /2TnV ■ (/5 2 u/T) . (33) 

The simplest evolution equations for the causal thermodynamic fluxes {77, q, n} in agreement 
with the second law (17) are again linear relationships, 

77 = -CX, (34) 

q = -ATY h , (35) 

7t = -2 V (Z)' h . (36) 
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8 J. Peitz and S. Appl 

Two additional thermodynamic coefficients Kk(p,n) had to be introduced in (31)-(33) as a 
consequence of the ambiguity involved in factoring terms which involve products ilq and 
7T • q in (28). Furthermore, (34)- (36) contain terms involving gradients of the and (3j. 
Since the K k are unknown a priori, these terms could be important even if the gradients 
themselves are small (Hiscock & Lindblom 1983). Finally, in (31)-(33) we neglected further 
contributions due to additional coupling terms between {H", q, n} and a, cu, which can be 
shown to exist in kinetic theory (Israel & Stewart 1979). 

The complexity of the full evolution equations (34)-(36) makes applications tractable 
only if certain simplifications are made. A particularly simple set of evolution equations 
results from the assumptions (Maartens 1997) 

= k = «i , (37) 
= «o = ol\ , (38) 
~ V • (ftu/T) , (39) 

where (37) reflects essentially the lack of knowledge on the while (38) neglects the cou- 
pling between heat flux and viscosity. Implications of (39) are multifold and need to be 
justified after a particular solution was found using a parametrization for (3j(p, n). The evo- 
lution equations resulting from (34)-(36) under the assumptions (37)-(39) are of covariant 
relativistic Maxwell-Cattaneo form, 

r (Dn) h + n = n, (40) 

n(Dq) h + q = q, (41) 
T 2 (Dn) h + n = 7T, (42) 

with the relaxation times Tj(p,n) given by 

ro = CPo , 7i = \T(5i , r 2 = 2 V (3 2 , , (43) 

and {II, q, jr} the (re-named) standard thermodynamic fluxes as in (23)-(25). In contrast 
to the algebraic constraint equations (23)-(25), the evolution equations (40)-(42) are first 
order partial differential equations, which assure that in the LRF the viscous bulk/shear 
stresses and the heat flux relax towards their standard limits {il, q, ft} on time-scales Tj. 
The relaxation times Tj follow in principle from kinetic theory, but can be estimated as mean 
collision times, 1/r ~ nEv, with E the collision cross section and v the mean particle speed. 



© 0000 RAS, MNRAS 000, 000-000 



3+1 formulation of non-ideal hydrodynamics 9 

For later use we re-write (40)-(42) as 



The conservation laws (5), (19), (20) and the Einstein equations (11), together with the 
evolution equations (34)- (36) constitute a complete system of hyperbolic first order PDEs 
for the solution vector of 24 (=1+1+3-1-10-1-1-1-3-1-5) dynamical variables {n, p, u, g, II, 
q, 7r}. This system represents a causal and stable theory for dissipative fluids (Hiscock & 
Lindblom 1983). 

2.6 Weakly dissipative fluids 

In many applications the inertia due to the dissipative contributions {II, q, n} can be ne- 
glected. In addition, it is convenient to simplify the evolution equations (40)-(42), which 
depend on the kinematic properties {0, a, cr} of the fluid, among which <y is particularly ex- 
pensive to compute. An appropriate simplification can be obtained by calculating (40)-(42) 
under the assumption of vanishing acceleration. Thermodynamic fluxes and thermodynamic 
forces calculated in this limit are underlined in the following. Note that the assumption of 
geodesic trajectories is only made for the calculation of the dissipative terms and not for 
the dynamics, i.e. we do not assume a geodesic velocity field satisfying Vuu = 0, but rather 
leave u unspecified. In the following this will be referred to as the weakly dissipative limit. 
In this limit the standard constraint equations (23)- (25) reduce to 




TO 
1 



(44) 




(45) 



(46) 



n = -(0, 



7T = — 2l]G , 



q = -ATDlnT, 



(47) 
(48) 
(49) 



with the thermodynamic forces calculated from the kinematic properties 



= 0\ a=o , a = a| a=0 = , a = a| a=0 = (Vu) s - l/30h . 



(50) 



The causal evolution equations (44)- (46) simplify to 




(51) 



(52) 
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10 J. Peitz and S. Appl 

Dn = -i (n - tt) . (53) 

3 DISSIPATIVE HYDRODYNAMICS IN 3+1 FORMULATION 

Particularly useful for time-dependent calculations in general relativity is the 3+1 formu- 
lation, where time derivatives are always with respect to globally defined universal time. 
Applications of 3+1 hydrodynamics and 3+1 magnetohydrodynamics have been mostly re- 
stricted to ideal fluids (see Bonazzola et al. 1993 for an exception). A collection of numerous 
general relativistic equations in the 3+1 representation can be found in Durrer & Straumann 
(1988). We give a 3+1 representation of relativistic dissipative hydrodynamics for both, stan- 
dard and causal thermodynamics. For a detailed derivation of the equations presented in 
this section we refer to Peitz (Peitz 1998). 



3.1 Generalities on the 3+1 formalism 

Assuming that spacetime (Ai,g) admits a slicing by slices S t , i.e. there is a diffeomorphism 
$:Ai^Ex/,7cE, such that the manifolds E t = $ _1 (E x {t}) are spacelike and the 
curves $ _1 (m, t) are timelike. These curves define a vector field d t which can be decomposed 
into normal and parallel components relative to the slicing, 

dt = ah + (3 . (54) 

Here n is the timelike unit normal field (congruence of fiducial observers=FIDOs) and (3 is 
tangent to the slices E t . a is the lapse function and |3 is the shift vector field. A coordinate 
system {x 1 } on E induces natural coordinates on M, i.e. $ _1 (m,t) has coordinates (t,x l ) 
if m G E has coordinates x l . The timelike curves 3 t have constant spatial coordinates 
(preferred timelike curves). Now set (3 = (3 l di (where 3j = d/dx 1 ). From g(n, 3;) = one 
finds g(3t, 3*) = —(a 2 — f3 l f3i) and g(3 t , 3j) = In coordinates co-moving with the FIDOs 
the metric thus reads 

g = -(a 2 - /JV^df ® dt + ftdf <g> dx* + (3 i dx i (g)dt + 7ydx i ®dx j (55) 
= -a 2 dt ® dt + 7y(d^ + (?dt) ® (dx j + (3 j dt) . (56) 

The forms dt and dx l + (3 l dt are thus orthogonal, y is the metric induced on E t , and the 
affine connection on (E t ,y) is denoted by V. 

The tangent and cotangent spaces of Ai have two natural decompositions, which give 
rise to two types of bases of vector fields and 1-forms. These are the dual pair {3 M } and 
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3+1 formulation of non-ideal hydrodynamics 11 

{da; M } for comoving coordinates {a+} and, on the other hand, the dual pair {rl, 5^} and 
{adt, dx 1 + f3 l dt}. Instead of the coordinate basis {3j} one may also use an orthonormal 
horizontal basis {e{\ with g(e i ,e :) ) = Sij, together with the dual basis {d 3 } for {dx*}. Then 
one has the two dual pairs {5 t ,e;}, {dt,£ 1 } and {e = n, e;}, {9 M } with the orthonormal 
tetrad {9^} given by 9° = adt and d i = & + /3*dt, with f3 i defined by |3 = (3%. From (54) 
follows the relation 

eo = n = - (3 t - P) . (57) 
a 

The 3+1 representation of respectively a scalar field A, a vector field B and a symmetric 
rank-2 tensor field C defined on (Ai,g) is understood by a representation with respect to 
the basis {eo = h, e^}, which we shall write in the following as 

A = E , (58) 

A 

B = Ee + S , S = S% , (59) 

B B B B 

C = £eo®e + S®eo + e ®S + T , S = S^i , T = T ij e>i <8> e,- . (60) 

c c cccccc 

The vector field Sb on (E t , y) corresponding to the vector field B on (M,g) is referred to as 
the horizontal part of B and, accordingly, the tensor field Tc on (^t> Y) is t ne horizontal part 
of C. Horizontal tensor fields appear bold face. They can be viewed as spatial components 
of the fields on (A4,g) (which appear as a subscript), after having been projected onto 
the 3-space orthogonal to h by ( ) y . If B = Sb, B is said to be a spatial vector field and, 
respectively, C is a spatial tensor field if C = To 

The 3+1 representation of the afline connection V on (M,g) depends on the connection 
forms, which depend on the kinematic properties of the FIDO congruence according to the 
irreducible decomposition of Vh in the sense of (4). The FIDO's kinematic properties are 
distinguished from kinematic properties of the fluid by a hat, i.e. 6", d>, a and are the 
shear, vorticity, acceleration and expansion of the FIDO congruence, respectively. The fields 
6", cu and a live in (E t , y) and are therfore spatial tensor fields, which we denote by 6" = T&, 
<i> = Ttb and a = Sa- The induced metric y on (£ t , y) is the horizontal part of the projector 
y = g + h ® h into the FIDO's frame. 

The FIDO world lines are orthogonal to the hypersurfaces and thus rotationsfree, uo = 
tu = 0. The FIDO's acceleration a is related to the lapse function by a = Vina. The 
connection forms on (Ai,g) can then be written in terms of the horizontal connection forms 
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12 J. Peitz and S. Appl 

on (E t , y) , and the horizontal parts of only two spatial tensor fields, namely 
a = Vina, (61) 
K = T • (62) 

K 

Here K is the extrinsic curvature tensor (second fundamental form), defined on (A4,g) by 

K = _vh = -i£ fl y. (63) 

An equation for K is obtained from the 3+1 representation of Einstein's equation (cf. (83) 
below). The 3+1 representation of the second equality of (63) is also recovered in the 3+1 
formulation of Einstein's equations (cf. (82) below), and provides an equation for y. This 
equation may be written in an alternative form based on the 3+1 representation of the first 
equality in (63), namely as 

K = -{v+\®y+^ t y) , (64) 

which confirms that K is indeed spatial (recall that d) = 0) and, furthermore, that K is 
symmetric, K = (K) s . 

For later use we give the 3+1 representation of DA, DB and DC according to (58)-(60), 

(65) 

DA 

E = - [d t E+VE-F- S , (66) 

DB 

S = - [d t -Cn)s+VS-EF- 1 K-S , (67) 

£>B G: V l/B BB B 



7 


(e,- 


a 


1 


(a,- 


a 


1 


(a,- 


a 




1 


(a ( - 


a 


1 


(«,- 


a 




1 


(a ( - 


a 



E = - [dt-Ca) E+VE-2F -S , (68) 
i)C <y \ ' / c c c 

S = - (dt-C^S+VS-lK-T- [ey + t] -F, (69) 

DC 

T = - (d t -Cn ) T+DT-2 7 K-T-F® S - S ®F , (70) 
DC a V l D J c c c c c 

with the horizontal vector field F defined on (S t ,y) by 

F = - 7 Vlna + K-u. (71) 

£p is the Lie derivative on (E t ,y) with respect to (3. The horizontal projection operator 
h = y + u ® u defined on (E t , y) projects into the 2-space orthogonal to u. Operators ( ) u 
and ( )h, T) = V u = u • V and D = Vh = h ■ V as well as ( ) s and ( ) a on (E t ,y) are 
to be understood in analogy to the corresponding operators on (Ai,g) ■ According to the 
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irreducible decomposition 



Vu = a+ cu + -i?h - a® u , (72) 

the kinematic properties of the fluid on (E t ,y) are given by tr = (Vu) s h , cu = (Vu)j^, 
a = Vu and i? = V • u. 



3.2 Conservation laws and Einstein's equations 

The velocity u has the 3+1 representation 

u = 7(eo + v) = 7e + u , (73) 

where 7 = E u = ~ n • u is the Lorentz factor with respect to the FIDOs and u = S u = 7v 
is the horizontal part of y ■ u = (g + h <g> h) • u = u + (h • u)h = u — 711. Since {e^} is 
orthonormal, the v l are physical 3-velocity components as measured by FIDOs. Similarly, 
the particle current n has the 3+1 representation 

n = n7(e + v) = n(7e + u) . (74) 
The stress energy tensor T in (15) has the 3+1 representation 

E = (p + P + i7)7 2 - (p + n) + 2 7 E + E , (75) 

T q 7t 

S = (p + p + i7)7u + 7S + £u + S , (76) 

T q q 7t 

T = (p + p + i7)u®u+(p + i7)y + u®S + S®u + T (77) 

T q q 7T 

= pu ® u + (p + i7)h + u®S + S®u + T • (78) 

q q 7t 

Note that q and n are orthogonal to u but not to h, and consequently q and n are no spatial 
fields. 

The 3+1 representation of particle number conservation (5), formally given by = E^. Q , 
can be written as 

= - (3 t -£ R V 7 n) + -V-(ami)-7ntr(K) . (79) 

The 3+1 representation of stress-energy conservation (6) splits into the energy equation, 
= -Ev-t, an d the momentum equation, = Sy-T, which can be regarded as the spatial part 
of = y ■ (V ■ T). Leaving T unspecified, these are 

= - f3 t -£ R ) J B-tr(K)E+2(S-V)lna + V-S-tr(K-T) , (80) 

CY \ I / T T T T T 

= - (d t -C R ) S + £Vlna-2K-S-tr(K)S+-V- (aT) . (81) 

a\l J /TT T t a T 

For completeness we give the 3+1 representation of the Einstein equations (e.g. Durrer 
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& Straumann 1988). They consist of evolution equations for y and K, obtained from the 
definition of the extrinsic curvature (63) and from the (space, space) components of (11). 
These dynamical equations are first order differential equations, 

^(dt-C^y = -2K, (82) 

-(dt-C^K = Ri(y) - 2K 2 + tr(K)K - 8tt T +47r(tr(K) - e)y - He(lna) , (83) 

where K 2 = K K, Ri( ) is the Ricci tensor and He( ) is the Hessian on (E t , y) , respectively. 
The (time, time) and (time, space) components of (11) yield the Hamiltonian and momentum 
constraints, 

R - tr(K) 2 - tr(K 2 ) = 16n E , (84) 

T 

V-K- Vtr(K) = 8ttS , (85) 

T 

where R is the curvature scalar on (E t ,y) . 
3.3 The standard constraint equations 

The 3+1 representation of the constraint equations (23)-(25) for the thermodynamic fluxes 
{il, q, 7t} depends on the 3+1 representation of the kinematic properties 0, a, a of the fluid, 
which can be calculated to 

E = W + #-7tr(K), (86) 
e 

E = 7 W + Thf-K u , (87) 

a 

S = 7 W + a, (88) 

a 

E = \{l 2 -l) (2^-tf + 7 tr(K)) + 7 (X> 7 - K u ) , (89) 
S = ^(7 2 -l)W + |(H/-2^ + 2 7 tr(K))u+i(D7 + 7 a) - (iK u y + K) -u 



2 

7 

f - 2 



1( 7 2_ 1)W+ i^u + D 7 + 7 a) - I(w + ^-7tr(K))u- (^K u y + k) • ^9p) 
(^W®u + u®w) -^W-7tr(K))h + <T-7K, (91) 



with time derivatives contained in 
W E 



i(a t -£p) 7 + (u.V)liia, (92) 
= ^ (a t -£p)u-2K-u + 7 Vlno; . (93) 
Decomposing DlnT finally yields the 3+1 representation of the constraint equations 



a 

W = 
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(23)-(25), 



E 
n 



-CE , 





E = -XT 

q 



s 

q 
E 



T 

ft 



-AT 



( 7 2 -l) ^ (d t - £p) InT + 7 X>lnT + E 
7 u^ (^-/^lnT + DlnT + S 



cr 

-2^S , 



-277 T 

cr 



(94) 

(95) 

(96) 
(97) 
(98) 
(99) 



3.4 The causal evolution equations 

For the sake of simplicity the following discussion will be restricted to the causal evolu- 
tion equations (44)-(46). Generalization to the full evolution equations (34)-(36) is straight- 
forward. The required 3+1 representations of D17, Dq, Dn follow readily from (65)-(70). 
Remains to decompose products q • a in (45) and n ■ a in (46), 

(100) 
(101) 



E = EE + S-S , 

q-a q a q a 

E = EE + S-S , 

7t-a 71 a 7T a 

S = ES + T-S ■ 

7t-a an tt a 



(102) 



The 3+1 evolution equations for the causal thermodynamic fluxes {En, E q , S q , E n , Sn, T n } 
then follow from (44)- (46) to 

(103) 



7 



a 



d t - £p) E+VE 



n n 



— ( E-E 

To v 77 n 



- (d t - Co) E+VE = —(E-E)+-fE+F-S, 
a V 1 / q q T'l^q <\ ' q- a q 

-(d t -C 6 )s+VS = -(s-s) + Eu + J EF + 7 K-S , 

CY V ' / q q 7~1 ^ q 1 ' q a 1 1 

— (d t — Cn) E+V E = -(. 



+2 7 £ +2F- S , 

T2 V ft 71/ 7T-a 7t 



(104) 
(105) 
(106) 



- (d t - £n)s+VS = — fs-s) + 7 S + £u + 7 K-S + ( Ey + t) F , (107) 

CtV ' / 7t 7t 72 ^ 1 11 ^ 71-51 7t " a 7T V 7t 7t / 

— (t-t)+u®S + S ®u + 2 7 K-T+F® S + S®F ,(108) 

To V ft 71/ 7t-a 7t-a 7T 7t 7T 



- (a t -£ 6 W+x>T = 

V l J / 7t 7C 

with F according to (71). 
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3.5 The weakly dissipative limit 

The limit of weak dissipation discussed in section 2.6 implies three simplifications of the 3+1 
equations of causal hydrodynamics. The first simplfication concerns the 3+1 conservation 
laws (80) and (81) for energy and momentum, where dissipative contributions {II, q, n} 
to the matter sources {Et, St, Tt} i n (75)-(78) can be dropped wherever they enter alge- 
braically, and only temporal and spatial gradients remain. The second simplification affects 
the 3+1 representation of the kinematic properties of the fluid, (86)-(91). Using the geodesic 
conditions E a = S a = 0, (87) and (88) yield simplified expressions for W and W, 

W = -X> In 7 + 1^/7 , (109) 
W = -a/7 , (110) 

which are no longer time-dependent. Therefore, any kinematic properties (86)-(91) simplify 
to time-independent expressions 

eE = $-V\n 1 + K u / 1 - 1 ti(K) , (111) 
& E = 0, (112) 

aS = 0, (113) 

S E = i(7 2 + 2)(i)ln 7 -^ u /7)-^(7 2 -l)(^-7tr(K)) 

= I (Z> In 7 - tf u / 7 ) + i (# - 7 tr(K)) + 1 (l> In 7 - K u / 7 - # + 7 tr(K)) 7 2 , (1 14) 

= i(a/ 7 + V 7 )+K-u + i(l)ln7-K u /7-^ + 7tr(K))7u, (115) 

aT = a-7K + i(x>ln7-K u /7 + 7tr(K))h, (116) 

with a = (Vu) s — l/2$h. Finally, the weak dissipation limit affects the 3+1 evolution equa- 
tions (103)-(108) for the thermodynamic fluxes {En, E q , S q , E n , Sn, T n }, where products 
Eq-s,, En-a, Sn-a as in (100)-(102) are dropped. This yields the weak dissipative evolution 
equations for the thermodynamic fluxes {En_, E q , S q , En, Sn, Tn}, 

-fa*- Cq) e+ve = - 

a V l J / 77 n r \ n n 

-fa*- Cq)e+ve = - 

a V V J q q Ti V q 



d t ~ E+VE = —[E-E], (117) 



9t-£ fi )£+X>£ = -f E-e) +F- S , (118) 

l/q q n \ q q / q 

- f3t-£ R ) S+I>S = -f S-s) +EF + 7 K-S , (119) 

ttV '*Vq q 7~i V q q / q q 

^ (d t - C a)E+V E = -( E-e) +2F-S , (120) 
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1 (^-Cq)s+VS = -fs-s) +7K-S + f Ey + t) F , (121) 

+ 2 7 K-T+F® S + S0F . (122) 

4 STATIONARY AND AXISYMMETRIC BACKGROUND SPACETIMES 

The 3+1 equations of dissipative hydrodynamics are specified to the class of stationary, ax- 
isymmetric background spacetimes. This situation is realized if the fluid under consideration 
has negligible influence on the gravitational field of a central object, and in addition this 
field is known to be stationary and axisymmetric. These assumptions include most appli- 
cations related to accretion/ejection flows in the vicinity of compact objects. For rotating 
black holes the vacuum metric is Kerr, and we give the equations also for this special case. 



2 

a V 



£a)T+VT = - 



1 

T"2 



T-T 

7t 7T 



4.1 Implications of symmetries 

The general form of a stationary, axisymmetric vacuum spacetime can be put in a form 
which is symmetric under a simultaneous change of sign of t and 0, the Killing coordinates 
associated with the commuting time and axial Killing vector fields k and m. Choosing the 
remaining two meridional coordinates as spherical coordinates allows to write g as 

g = - a 2 dt <g> dt + lu 2 (d0 - Lddt) <g> (d0 - udt) + e 2 ^dr <g> dr + e 2v A6 <g> d# , (123) 

with the invariant metric coefficients 

Co 2 = m 2 , w = -k-m/m 2 , a 2 = -k 2 + k • m/m 2 . (124) 

For the physical interpretation of u, u and a see e.g. Bardeen (1970). The generic choice of 
the fiducial congruence is (see e.g. Thorne & Macdonald (1982) for criteria that uniquely fix 
this choice) 

h = eo = -(k + um) . (125) 

a 

These FIDOs possess vanishing specific angular momentum, = h • m = 0. Therefore they 
correspond to Bardeen's (1970) zero angular momentum observers (ZAMOs). Furthermore, 
since h • m = 0, m is a spatial vector, which we denote on (E t , y) by m = S m - Comparison 
with (57) shows that the shift vector has to be chosen as 

(3 = -cum , (126) 
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and the metric induced on (S t ,y) becomes 

y = u 2 d(j) <g> d(j) + e 2 "dr <g> dr + e 2iy d# <g> d6 . (127) 

The 3+1 representation of Killing's equation for k and m, together with their commutivity, 
allows to establish the following relations (Thorne & Macdonald 1982) 

d tf i = , d t u = , d t u = , d t m = , d t y = 0, (128) 
m-V/i = 0, m-Vz/ = 0, m • Vcu = , (Vm) s = , (129) 
= 0, cr = -J- ( m <g> Vw + Vw <g> m) . (130) 

Note that (Vm) s = states that m is a Killing vector field on (S*,y) . However, as a 
consequence of 'Vuj ^ 0, |3 is not a Killing vector field on (S 4 ,y) . As a concequence of 3 t y 
and = 0, the horizontal field K according to (64) reduces to 

K = -a = --J-fm® Vu + Vcu^m) . (131) 
Therefore K measures the shear of hypersurfaces S t , which vanishes for = 0. 



4.2 Gravitomagnetic and gravitoelectric tensor fields 

A characteristic phenomenon in axisymmetric spacetimes is the dragging of inertial frames. 
Physical implications due to this effect (see Thorne et al. 1986 for the case of the Kerr metric) 
are often described in terms of the gravitomagnetic tensor field H, defined on (S t ,y) by 

H = -V(3 = --V(wm) = -- [wVm + Vuj® m] . (132) 

a a a 

It is clear from (131) that K = (H) s , and therefore 

K-T = - (H • T + T • H) . (133) 

The antisymmetric part of H can be expressed as an axial vector field on (S t ,y) , namely 
the gravitomagnetic vector field 

J = -VA0. (134) 
a 

For an arbitrary vector field S on (S t ,y) one finds the relation 

H-S = K-S-ijAS. (135) 

The Lie derivatives £p in the 3+1 equations derived in section 3 can then be expressed as 

C$E = (|3-V)£, (136) 
£pS = (0-V)S-(S-V)0 = (0-V)S-aS-H, (137) 
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£pT = (0- V)T-2((T- V)|3) s = (0 • V)T - 2a(T • H) s = (0 • V)T - 2T • K . (138) 

The term S • H in (137) may be expressed by either H as in (132) or in terms of J and 
K via (135). In addition to the gravitomagnetic tensor fields H and J, we introduce the 
gravitoelectric vector field G, defined on (E t ,y) by 

G = -a = -Vina. (139) 
This field measures the gravitational acceleration measured by the fiducial observers. 

4.3 Conservation laws 

The 3+1 conservation laws for particle number (79), energy (80) and momentum (81) in a 
stationary, axisymmetric background read 

= - ( |- - • vV 7 n) + - V(cmu) , (140) 
a \ot J a 

= - (!-- 0- V ) S+V-S-2S-G- tr(K-T) , (141) 

a \Ot J T T T T 

= - 0- v) S-£G-H-S+-V(aT) 

a \ot J t t t a t 

= -f^-P-V^S+V-T-H-S-fsY + TVc. (142) 

a\Ot / T T T \T T/ 

4.4 Evolution equations 

The 3+1 representation (86)-(91) of the kinematic properties of the fluid in stationary, 
axisymmetric background read 

E = W + $, (143) 
e 

E = 7 W + i> 7 + 6- u , (144) 



S = 7 W + a, (145) 
( 7 2 -1)(W- + 7 (X> 7 + o- u ) , (146) 



a 

3 

f = 2 



5 = i(7 2 -l)W+^(^-2^u + i(D 7 + 7 a)-(6- u y+6-)-u, (147) 

T = J(w<g>u + u<g> w) - ^WTi+ cr + 7 6- (148) 

6 2 V / 3 

with time derivatives contained in 

W = IfA-p.vV-^u, (149) 



ct \3t 

W = -(|--0Vlu-Hu- 7 G. (150) 
a \ot 
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The 3+1 representation of the standard constraint equations (94)-(99) can be written as 

(151) 

(152) 



E 
n 



<E , 
e 



E = -kT 

q 



-kT 



1 /3 



( 7 2 -l) -[jj- • V ) bl ~ VluT - / 
7U- ( |- - |3 • V I InT + DlnT + S 

a \Ot I a 



S = 

q 

E = -%VE , 

ft cr 

S = -2r7S , 

ft cr 

T = -2r]T ■ 

ft cr 

The 3+1 representation of the causal evolution equations (103)-(108) reads 

l(l 
a\dt 



(3 • V E +T> E 

n n 



(e-e) , 

V fi n / 



E 

TO V 77 



-(|r-P-V)S+I>S = — (s-s) + £u + £F + ^7JAS, 

CH \Ot J q q T"i V q q / q-a q 2 q 



if! 

l(l 
a \dt 

1(1 

a \dt 

1(1 

a \dt 

with 



P • V ] E+VE 

q q 



Ti V q q 



— [E-E)+lE +F-S , 



q-a 



(3 • V )E+VE 

J Tt Tt 

- • v) s+x>s 

J Tt Tt 

• V ) T +t> T 



— f + 2 7 £ +2 7 F- S , 



S-s)+7S+£u + i 7 JAS+(£Y + T) • F ,(161) 

7~2 V ft 7t/ 7T-a 7t-a 2 



7t 7t 



!(' 



Tt \ 7t 



(153) 
(154) 
(155) 
(156) 

(157) 
(158) 
(159) 
(160) 



T-t) +u® S + S ®u + 7 F® S+7S®F , (162) 

T2 V ft Tt J 



7t-a 7t-a 



71 7t 



7 G — (X • U . 



(163) 



4.5 Weakly dissipative limit 

The 3+1 representation of the kinematic properties of the fluid in the weak dissipation limit 
can be written as 



e_E = d - £>m 7 + 6r u / 7 , 
*E = 0, 
aS = 0, 



oE = i( 7 2 + 2) (l)ln 7 -6- u / 7 ) -i( 7 2 -l)^ 

(z>ln7-a- u / 7 ) + ^ + i(x)ln 7 -or u / 7 -^ 7 2 , 



2 
3 



(164) 
(165) 
(166) 

(167) 
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2 S = ^(a/7 + V7) + &-u + ^(pin 7 -CT u /7- 0)711, (168) 

aT = a- 7 a+^Plii7-CT u /7)h. (169) 

The 3+1 weak dissipation evolution equations (103)-(108) for {En_, E q , S q , En, Sn, Tn} can 
be written as 

1 (^7-P-v)e+ve = -(e-e), (170) 

a \ot / n n r \ fr R ' 

1 (^7- P-V) E+VE = -( E-E) +F-S , (171) 
a \ot /qqT-iVqq/ q 

- (It- P-V)S+I>S = —( S-s) +EF + ^3 AS , (172) 

a \0t J q q Ti\qq/q Z q 

-(^7-&-v \ E+VE = —( E-E) +2F-S , (173) 

Q: \ OT y n n T2 \ n n / n 

~ (Il- I 3 ' V ) S+X>S = — fs-s) + ^7JAS + f J By + T N ) F, (174) 
a yot y n 7t T2 \ ft n / 2 ^ z* ' 

- (It- P- v ) t+x>t = — (t-t) +f® s + s®f (175) 

with F according to (163). 

4.6 Specification to Kerr geometry 

The behaviour of fluids in the vicinity of a rotating black hole is governed by the conservation 
laws and evolution equations in Kerr metric. The Kerr metric is a two parameter family of 
stationary, axisymmetric vacuum spacetimes. Parameters are the mass M and the specific 
angular momentum a = J/M of the black hole. In terms of the functions 

A = r 2 + a 2 -2Mr, q 2 = r 2 + a 2 cos 2 9 , E 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 6 (176) 

the lapse function a and the non- vanishing components of the shift vector (3 and the metric 
y according to (56) are 



~2 



a = — , ^ = - UJ = -^ r , ^ = ^sin£ , e 2 ^ = ^ , e 2v = q 2 . (177) 

A consistent treatment of viscous hydrodynamics in the vicinity of a rotating black hole 
would require a careful analysis of the boundary conditions at the horizon. Such an analysis 
is best performed within the concept of a streched horizon, defined by a small value of lapse 
a. We refer to the membrane paradigma of Thorne et al. (1986), where this approach is 
applied to ideal magnetohydrodynamics and Maxwell's equations. This work also contains 
a coordinate representation of G and H. 
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4.7 Stationary and axisymmetric flows 

The assumptions of stationarity and axisymmetry so far concerned exclusively the back 
ground metric. In many problems it is justified to assume these symmetries to hold also for 
the fluid configuration under consideration. The equations for stationary flow are obtained 
by setting d t = in the 3+1 equations, since now k is a Killing vector field also for the fluid. 
Similarly, the equations for axisymmetric flow follow from setting (3 • V = 0, since m then 
is a Killing vector field also or the fluid. 

The post-Newtonian limit of causal viscous hydrodynamics is recovered from the 3+1 
equations by setting a — 1 (i— > G = 0) and (3 = (i— > H = a — 0), and furthermore setting 
7 equal to unity, while gradients of 7 are retained. This limit contains still non-vanishing 
quantities of post-Newtonian order V u u. Neglecting these then yields the Newtonian limit 
of causal viscous hydrodynamics, where Ee, En, E a , So- vanish. A related causal description 
of viscous angular momentum transfer in Newtonian accretion disc boundary layers was 
considered by Papaloizou & Szuskiewicz (1994) and Kley & Papaloizou (1997). 

5 CONCLUSIONS 

We have provided a complete set of equations for dissipative relativistic hydrodynamics in 
their 3+1 representation. Furthermore, we have specified the general system to the class 
of stationary axisymmetric vacuum spacetimes, with the Kerr metric as the most relevant 
astrophysical example. For this case we have written the equations in a form where the 
dragging of inertial frames is described by the gravitomagnetic tensor field. This allows to 
combine the equations with the 3+1 Maxwell's equations as in Thorne & Macdonald (1982). 

Causality has been accounted for by using the extended causal description of thermody- 
namics, relativistically formulated on a phenomenological level by Israel (1976) and Stewart 
(1977) and justified by kinetic theory (Israel & Stewart 1979). In contrast to the convention- 
ally used compressible Navier-Stokes description of non-ideal hydrodynamics, the equations 
of extended causal thermodynamics guarantee finite propagation speeds of heat and viscous 
signals and yield stable local thermodynamic equilibria (Hiscock & Lindblom 1983). 

A causality preserving formulation is required whenever the thermodynamic timescale 
becomes comparable to the dynamical timescale and therefore the assumption of local ther- 
modynamic equilibrium is not justified. This is particularly the case in supersonic flows 
and/or processes in the vicinity of the event horizon. Astrophysical examples and thus po- 
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tential fields for application of causal thermodynamics in the formulation presented here 
include the gravitational collapse of stars (e.g. Baumgarte et al. 1995), the innermost parts 
of accretion discs around black holes (e.g. Peitz & Appl 1997) or the collision of neutron 
stars (e.g. Rasio & Shapiro 1992). In many problems of interest the inertia due to the dis- 
sipative contributions to the stress-energy tensor can be neglected. The 3+1 formulation of 
the corresponding simplified set of equations is given. 

The five conservation laws for particle number, energy and momentum, the ten Einstein 
equations for the metric tensor and the ten evolution equations for the thermodynamic 
fluxes form a hyperbolic system of first order PDEs tractable by numerical methods (e.g. 
Bonazzola et al. (1993)). 

We are grateful to M. Camenzind, R. Khanna and S. Spindeldreher for enlightening 
discussions. Support by Deutsche Forschungsgemeinschaft under SFB 328 (JP) and SFB 
359 (SA), and by the French Ministery of Foreign Affairs (SA) is acknowledged. 

REFERENCES 

Bardeen J.M., 1970, ApJ, 162, 71 

Baumgarte T.W., Shapiro S.L., Teukolsky S.A., 1995, ApJ, 443, 717 
Bonazzola S., Gourgoulhon E., Salgado M., Marck J.A., 1993, A&A, 278, 421 
Durrer R., Straumann N., 1988, Helv. Phys. Acta., 61, 1027 
Eckart C, 1940, Phys. Rev., 58, 919 

Gammic C.F., Popham R., 1997, submitted to ApJ (astro-ph/9705117) 

Hiscock W.A., Lindblom L., 1983, Ann. Phys., 151, 466 

Israel W., 1976, Ann. Phys., 100, 310 

Israel W., Stewart J.M., 1979, Ann. Phys., 118, 341 

Jou D., Casas-Vazquez J., Lebon G, 1997, Extended irreversible thermodynamics. Springer Verlag, Berlin 
Kley W., Papaloizou J.C.B., 1997, MNRAS, 285, 239 

Landau L.D., Lifshitz E.M., 1959, Fluid Mechanics. Pcrgamon Press, London 
Levermore CD., Pomraning G.C., 1981, ApJ, 248, 321 
Lindblom L., Hiscock W.A., 1983, ApJ, 267, 384 

Maartens R., 1997, in Maharaj S.D., ed., Proceedings of Hanno Rund Workshop. Natal Univ. Press, to appear (astro- 
ph/9609119) 

Maartens R., Triginer J., 1997, submitted to Phys. Rev. D (gr-qc 9707018) 
Macdonald D., Thornc K.S., 1982, MNRAS, 198, 345 
Miiller I., 1967, Z. Physik, 198, 329 
Narayan R., 1992, ApJ, 394, 261 

Papaloizou J., Szuskicwicz E., 1994, MNRAS, 268, 29 

Peitz J., Appl S., 1997, MNRAS, 286, 681 

Peitz J., 1998, PhD thesis, University of Heidelberg 



© 0000 RAS, MNRAS 000, 000-000 



24 J. Peitz and S. Appl 



Popham R., Narayan R., 1992, ApJ, 394, 255 

Rasio F. A., Shapiro S.L., 1992, ApJ, 401, 226 

Stewart J.M., 1977, Proc. Roy. Soc. (London) A., 365, 43 

Syer D., Narayan R., 1993, MNRAS, 262, 749 

Thorne K.S., Macdonald DA., 1982, MNRAS, 198, 339 

Thorne K.S., Price R.H., Macdonald D.A., Sucn W., Zhang X., 1996, in Thorne K.S., Price R.H., Macdonald D.A., eds., Black 

Holes: The Membrane Paradigm. Yale Univ. Press 
Zimdahl W., Pavon D., Maartens R., 1997, Phys. Rev. D, 55, 4681 (aph-ph/9611147) 



© 0000 RAS, MNRAS 000, 000-000 



